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Irreducible Many-Body Casimir Energies of Intersecting Objects 
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PACS 1 1 . 10 . -z - Field theory 
PACS 1 1 . 10 . Gh - Renormalization 

PACS 42.50.Lc - Quantum fluctuations, quantum noise, and quantum jumps 
PACS 1 1 . 80 . La - Multiple scattering 

Abstract. - The vacuum energy of a bosonic field interacting locally with objects is decomposed 
into irreducible iV-body parts. The irreducible TV-body contribution to the vacuum energy is finite 
if the common intersection 0\ n O2 ■ ■ ■ n On of all N objects Oi, i = 1, . . . , N is empty. I prove 
that the perturbative expansion of the corresponding irreducible A-body spectral function <f>' N ' (/3) 
for P ~ vanishes to all orders even if some of the objects intersect. These irreducible spectral 
functions and their associated Casimir energies in principle can be computed numerically or ap- 
proximated semiclassically without regularization or implicit knowledge of the spectrum. They 
are analytic in the parameters describing the relative orientation and position of the individual 
objects and remain finite when some, but not all, of the N objects overlap. The Feynman-Kac 
theorem is used to compute Casimir energies of a massless scalar field with potential scattering 
and the finiteness of iV-body Casimir energies is shown explicitly in this case. The irreducible 
iV-body contributions to the vacuum energy of a massless scalar field with potential interactions 
is shown to be negative for an even- and positive for an odd- number of objects. Some simple 
examples are used to illustrate the analyticity of the iV-body Casimir energy and its sign. A mul- 
tiple scattering representation of the irreducible three-body Casimir energy is given. It remains 
finite when any two of the three objects overlap. 
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Introduction. — The Casimir energy for two disjoint 
bodies is finite and may be estimated [IrE] . It can, in prin- 
ciple, be computed to arbitrary numerical precision [BHH]. 
For disjoint bodies, the multiple scattering representation 
of the interaction energy [9HI2] thus solves many problems 
encountered in technological applications [6l[T3]. We here 
develop an extension of this formalism and extract finite 
irreducible Casimir energies for more than two bodies that 
are not necessarily mutually disjoint. The analysis gives 
a new interpretation to finite parts of zero-point energies 
that could provide a framework for exploring gravitational 
effects due to vacuum energies [14j and result in a more 
systematic approach to Casimir self-stresses for arbitrarily 
shaped bodies. 

For clarity of presentation and to avoid infrared issues, 
we assume that the objects {Oi\ i = 1, . . . , N} are all em- 
bedded in a large, but finite, connected Euclidean region 
2)0 of dimension d. The thermodynamic limit 2>0 — > lZ d 
may be taken at the end. Formally, the vacuum energy 
£12. ..n due to a massless bosonic field in the presence of 



N objects can be decomposed into 

£ 1 2...n = ^+J2^ 1) +J2^ 



?(2) 
ij 



i<j 



-£ {N) (1) 



where £^ , is the irreducible contribution to the vacuum 
energy that depends on all k objects . . . Oi k in the do- 
main 2)0. Eq. (fTJ) would recursively define the irreducible 
iV-body Casimir energy N as an alternating sum of 
vacuum energies (see eq. (fT12j) ). The irreducible two-body 
Casimir energy for instance is 



?(2) 
'12 



£1 — £2 + £$ 



(2) 



It is finite for two disjoint objects. Below I show that 
the irreducible TV-body Casimir energies £[2^ n are S^te 
as long as the common intersection 0\V\Oi ■ ■ -PiOn of all N 
objects is empty. For two objects this requires them to be 
disjoint, but three and more objects need not be mutually 
disjoint and the irreducible three-body contribution to the 
Casimir energy of a triangle, for instance, is finite. For a 
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massless scalar field whose interaction with the objects is 
modeled by positive local potentials, the irreducible TV- 
body Casimir energies are shown to be negative for an 
even-, and positive for an odd-, number of objects. 

Subtracted TV-body Spectral Functions. Some 
irreducible vacuum energies, such as self-energies may 
diverge when the associated objects overlap, eq. (JI]) thus is 
formal in the sense that it requires a high-frequency cutoff. 
However, spectral functions generally are finite and well 
defined even when the one-loop vacuum energies are not. 
They can be similarly decomposed into irreducible parts 
and we therefore relate the irreducible TV-body Casimir en- 
ergy £( N > to the corresponding irreducible TV-body spec- 
tral function (Ar) (/3), 

4>^ N \(3) is constructed as follows. Let 55 s represent the 
domain 23 with objects {Oj\j G s} embedded, T>i...n 
being the finite domain 55g with all TV objects included. 
Denote with SP(s) the power set of the elements of a set s 
of finite cardinality |s| < TV with <# N = qj({l . . . TV}). Let 
4> s (j3) be the spectral function, or trace of the heat kernel 
M$) a , for the domain 55 s , 

MP) = TrJW/3) = e- px »W*. (4) 

n(EN 

Here {An(55 s ) > 0, n € IN is the spectrum of a bosonic 
field that vanishes on the boundary of 550 and whose in- 
teractions with the objects in 55 s are local. We assume the 
interaction of the field with the objects may be described 
by positive local potentials or take the form of (compati- 
ble) local boundary conditions. 

The irreducible spectral function (f>( N ) (/3) of eq. (J3]) is 
the alternating sum of spectral functions </> s (/3) for the 
individual domains D Sy 

W (/*):= E (-1)" H,| ^G8) • (5) 

This is a special case of the geometrical subtraction pro- 
cedure advocated in ref. [15] . We will see that in this case 
the asymptotic expansion of the heat kernel vanishes to 
all orders if the common intersection of all TV objects is 
empty. A pictorial representation of eq. ([5]) for four line 
segments as objects in a bounded 2-dimensional Euclidean 
space is given in fig. [TJ 

To facilitate proving that the integral in eq. ([3]) is finite, 
we demand that the individual heat kernels are uniformly 
bounded by the free heat kernel of lZ d , 

< £ Ss (x, y; (3) < K(2ir/3)- d / 2 e- {x - y)2/m , (6) 

for some finite K > 0. For a scalar interacting with lo- 
cal positive potentials this is implied by the Feynman-Kac 





Fig. 1: The subtracted spectral function 4>- N '{j3) defined in 
eq. © for a bounded two-dimensional domain with four 
intersecting line segments as objects. Each pictograph repre- 
sents the spectral function of the corresponding domain taken 
with the indicated sign. Various local features that contribute 
to the asymptotic expansion of each spectral function at high 
temperatures (small f$) have been highlighted: lines of different 
color correspond to possibly different, but compatible, bound- 
ary conditions or local potentials. Since the intersections of 
line segments generally differ, each vertex is shown in a differ- 
ent color. The contribution to the asymptotic expansion from 
any particular local feature vanishes: the total signed num- 
ber of times any particular line segment contributes is zero, as 
is the total signed number of times any particular vertex oc- 
curs. A random path that crosses three of the four segments 
is shown schematically. Its contribution to the Casimir energy 
of a scalar field in the Feynman-Kac path integral vanishes. 

theorem |16) : the heat kernel is just the transition proba- 
bility for Brownian motion and a (positive) potential re- 
duces it. Dirichlet boundary conditions in particular, may 
be imposed on a surface by killing any random path that 
crosses it. Because any condition is satisfied only by a 
reduced number of them, the bound also should hold for 
objects represented by any other local boundary condition. 
For the following it is sufficient that correlation functions 
vanish faster than any power of /3 as /? — > for any finite 
separation (x — y) 2 > S 2 > 0. One thus may be able to 
relax the uniform bound of eq. © considerably. 

We assume that the spectrum is discrete and positive 
semi-definite. 4> S {I3) may be interpreted as a bosonic single 
particle partition function at inverse temperature (3 and a 
positive spectrum is equivalent to the absence of tachyons 
in the causal local theory. The spectral functions <j> s {0) 
of eq. (J5J) in this case are positive and monotonically de- 
creasing, approaching, at most, a finite positive constant 
for /3 ~ 00. Although we only treat a scalar bosonic sys- 
tem, the following also holds for electromagnetic fields in 
the absence of free charges. 

In local field theories, the asymptotic expansion of <j) s (/3) 
for small f3 has the general form [TTH^U] . 

00 

0,(0 ~ 0) ~ J2 Wr 72 ^ + 0(e- £ ™/< 2 «) , (7) 

u— — d 

where the Hadamard-Minakshisundaram-DeWitt-Seeley 
coefficients for the domain 55 s have length-dimension 
(— v). Note that if eq. (0 holds, exponentially suppressed 
terms are associated with classical periodic paths of finite 
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length t m i n . We decompose the heat kernel coefficients 
of eq. ([7]) into parts arising from local features of the 
individual objects and their overlaps, 

4"' = E ^ > (8) 
t£<P(s) 

where the sum extends over all (\s\\) sets in the power 
set ^3(s) of the set s. eq. © recursively defines reduced 

heat kernel coefficients ai"^: the a^' are the heat ker- 
nel coefficients associated with the Euclidean domain D^; 
the a/X give their change when object j is inserted; the 

{jk} accoun t f° r further changes in the asymptotic heat 
kernel coefficients due to local overlaps of objects j and 
k. Note that the are not the heat kernel coefficients 
of the domain D T — they are their irreducible part only 
and arise from arbitrary short correlations near common 
intersections of the objects in the set r. a^^y = for 
two disjoint objects j and k, if we assume (as implied by 
eq. ©) that asymptotic correlations over finite distances 
|x — y| > 6 > vanish faster than any power in ft. Simi- 
larly, £1/^23} = if 0\ n O2 n O3 = 0. = vanishes 
even if the three objects are not mutually pairwise dis- 
joint, a^jfr-j accounting for contributions to the asymptotic 
power series of the pairwise intersection Oj H Ok ■ 

This argument may be extended to N objects to imply 
that for local interactions the correction, 

a ti...N} = ' if O^-'-DOn =0 . (9) 

Note again that the condition in eq. © does not imply 
that the objects have to be mutually disjoint (except for 
N = 2). It then is a combinatoric matter to prove that 
for t C {12 . . . N} the contribution of any non-zero to 
the asymptotic expansion of <f>( N } (ft) in eq. ([5]) vanishes. 
Because the other \s\ — |t| objects may be selected from the 
remaining N— |r| in any order, the number of times the set 
r occurs as a subset of the sets in tyjy (with cardinality 
|s| > |-7-|) is the combination ^J^y.^l^y. = 
For TV > \t\ the contribution to the asymptotic expansion 
of (ft) in eq. ([5]) proportional to then is 

<2nfJ)»/*aP jr (-1)^1^-1^=0. (10) 

|.|=|r| VV |S|y 

When N objects have no common intersection, the asymp- 
totic expansion of the irreducible ./V-body spectral func- 
tion 4>( N } (ft) thus has the form, 

^»(^O)~0(e-^/( 2 «), (11) 

and vanishes faster than any power of ft. This may be ex- 
plicitly verified in examples like that shown in fig.[U noting 
that contributions to the asymptotic expansion propor- 
tional to the volume, surfaces, corners, curvatures etc. . . , 



all cancel. Together with the fact that the spectral func- 
tions 4> s (ft) decay monotonically and remain bounded for 
large ft, the asymptotic behavior of eq. (fTTj) implies that 
the Casimir energy given by the integral in eq. © is finite. 

The subtraction procedure allows one to formally inter- 
pret £~( N ) as the alternating sum of vacuum energies £ s 
associated with the domains D s , 

£W = £ (-1)^-1^ . (12) 

The sum on the right side of eq. (1121) requires some regu- 
larization to be meaningful but, if this procedure does not 
explicitly depend on the specific domain 2) s (for instance 
by regularizing the proper time integrals), the previous 
considerations show that the irreducible iV-body contri- 
bution £ ( N > remains well defined as the regularization is 
removed. The absence of a power series in the asymp- 
totic expansion of <j^ N \ft ~ 0) also explains why a semi- 
classical approach based on classical periodic orbits tends 
to approximate Casimir energies fairly well [TH3"l l2"TH2"3"] : it 
reproduces the leading exponentially suppressed terms of 
the asymptotic expansion. 

Massless Scalar Field with Local Potential In- 
teractions. — The subtraction procedure we have just 
outlined is particularly transparent for a massless scalar 
field in a bounded Euclidean space D® whose interaction 
with the objects is described by a local (positive) poten- 
tial V = Y2i Vi- It allows for an alternative proof in this 
case and provides an additional insight into the sign of 
the irreducible ./V-body contribution. Using the world-line 
approach of [7] for potential scattering, the Feynman-Kac 
theorem [16] generally states that, 

r d x 

*'W = h t wFF v »-toW> (13) 

where 'Pa, [j?3(x)] is the probability for a standard Brow- 
nian bridge], ^a(x) = {x t ,0 < t < /3;x = x t = x}, that 
starts at x and returns to x after "time" ft, to not exit 
35 s and survive its encounters with the objects. The sur- 
vival probability of any particular Brownian bridge in 35 s 
is given by p a (£p(x)) = exp[- /„ V s (x t )dt], where V s is 
the sum of local potentials representing the objects in 
Dirichlet boundary conditions in particular may be im- 
posed by setting p s — for a loop that crosses the surface 
of the object and p s — 1 if it does not. 

The contribution to 4> {N) (ft) of a loop ^ t) (x) that 
remains within iDg and encounters all objects of r C 



-"-A standard Brownian bridge £p(x) = {x + y/]3(W(t) — 
f.W(l));0 < t < 1} is generated by a standard oi-dimcnsional 
Wiener process with stationary and independent increments for 
which W(t > 0) is normally distributed with variance td and van- 
ishing average. 
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{1 . . . N} and no others is, 

i&V{t) s=| 7 | v 1 " 7 

s(1t = 7 

(14) 

The contribution vanishes independent of the survival 
probabilities p 1 . Only loops that touch all N objects, 
(r = {1,...,TV}) contribute to the alternating sum in 
eq. (J5J) and we have that 

W (« = (-1) JV / T^b^W], (15) 

where 'P'- Ar - ) [x; /3] is the probability that a standard Brow- 
nian bridge starting at x and returning to x after "time" 
ft does not exit 2)0 and is killed by every one of the N 
objects. Composing the probabilities one observes that a 
Brownian bridge is killed by every one of N objects with 
probability, 

^(killed by every one of N objects) = (— l)' 7 'p 7 . 

76'Piv 

(16) 

cq. (IT6l) is the extension to N objects of the statement, 

^(killed by 0\ and killed by O2) = 
= (1 - pi) + (1 - pa) - (1 - P12) =P<»-Pi-P2+ Pi2- 

Note that Eqs. (|15[) and (p~6|) do not require survival prob- 
abilities to be independent, P12 = p\Pi ~ which is the case 
only for potentials that do not overlap (and therefore cor- 
responds to mutually disjoint objects). Since is a 
positive probability, the factor of (— 1) N in eq. (fT5|) deter- 
mines the sign of 4>^ N \f3). [For Dirichlet boundary con- 
ditions on the objects, it arises because paths that touch 
all N objects contribute only to 4>@(13)-] For scalar fields 
interacting by local potentials with N objects that have 
no common intersection, the irreducible TV-body Casimir 
energy is finite and its sign is given by, 

(-l)^W < 0. (17) 

It is remarkable that the sign of £^ N > depends only on the 
number of objects and the fact that the interaction with 
the scalar field is described by local potentials. The irre- 
ducible scalar two-body Casimir energy in particular, in 
this case is negative independent of mirror symmetry |11) . 
This also holds for Dirichlet boundary conditions. But 
eq. (JTTJ) does not hold for boundary conditions like Neu- 
mann's, that are not described by potentials [30] ■ Also, 
the sign of the TV-body Casimir energy does not of itself 
determine whether Casimir forces are attractive or repul- 
sive. The subtraction clearly exhibits the irreducible part 
of the vacuum energy computed in eq. (fT2|) . but it is impor- 
tant to correctly interpret this energy. The finite TV-body 
Casimir energy obtained here is the irreducible TV-body 
correction to the vacuum energy that remains when all 



M-body vacuum energies with < M < TV have been 
accounted for. The latter may themselves be finite but 
very often are not, and the sign of £^ N ^ given in eq. (JTTJ) 
is that of the irreducible TV-body part only, which, in gen- 
eral, does not coincide with the work needed to assemble 
the TV objects. 

For a scalar field, eq. (fT5)l interprets 4>^ N \f3) as a prob- 
ability for a random walk that satisfies certain geometric 
conditions. Since they have to touch TV objects that have 
no common intersection, Brownian bridges that contribute 
in eq. (|15p necessarily are of finite length. The probabil- 
ity [x; /?] thus is bounded from above by the shortest 
closed classical path of length £ m [ n that just touches all 
objects, 

0<pW[x;/3]<e-^"/( 2 « . (18) 

For a domain of finite volume, the bound of eq. ([15} 
implies that the asymptotic power series in j3 of (fj ~ 
0) vanishes to all orders in this case, as we have previously 
argued more generally. 

Examples. — Consider the example of a scalar field 
in lZ d satisfying Dirichlet boundary conditions on (d + 1) 
intersecting, (d— l)-dimensional hyper-planes. In this case 
£ ( d + 1 ) indeed is the work required to adiabatically move 
the last hyperplane into position from infinity: van- 
ishes as the volume enclosed by the hyper-planes becomes 
infinite and depends continuously on their position. These 
are simple consequences of the smoothness and continu- 
ity of the probability that Brownian bridges cross all of 
them in time /?. eq. (|12p implies that (infinite) hyper- 
planes forming a simplex, such as a triangle(d = 2) or a 
pyramid(<i = 3), tend to repel (triangle) in even and to at- 
tract (pyramid) in odd dimensional spaces. Contributions 
to Casimir energies from interior modes of domains with 
generalized reflection symmetries have been computed an- 
alytically as well as numerically [24-26! . None of these 
obtain only the finite irreducible part of the three-body 
Casimir energy and the results are somewhat ambigu- 
ous because corner divergences are subtracted in differ- 
ent ways. The world-line method [7 HTM5 1] outlined above 
could provide fairly accurate numerical estimates for the 
irreducible TV-body part of the Casimir energy of a scalar 
field with Dirichlet boundary conditions on any generic 
set of TV intersecting hyper-planes without restricting to 
contributions from interior modes. The irreducible three- 
body Casimir energies of some weakly interacting inter- 
secting objects are explicitly found to be positive and fi- 
nite [35] . 

The analytically tractable Casimir energy of 2 d pairwisc 
parallel (d— l)-dimensional hyper-planes forming a multi- 
dimensional tic-tac-toe-like pattern in lZ d that encloses 
an inner hyper-rectangle with dimensions £1 X • • • X l& was 
previously considered in [27) . without separation of the 
analytic irreducible 2 d -body contribution. The irreducible 
2 d -body Casimir energy £ r 2 0Ct for a scalar field satisfying 
Dirichlet conditions on all the hyper-planes in fact has the 
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simple form 



hcY[{d+l)/2] 

47r (d+l)/2 



E-E 



(19) 



where V rcc t. = FIj=i IS the volume of the hyper- 
rectangle and L(n) = \jYl!j=i ^ s half the length of 
a classical periodic orbit in its interior that reflects up- 
times off the j-th pair of parallel hyper-planes. Only 
classical periodic orbits that touch all hyper-planes con- 
tribute to £j eci ■ Consistent with the previous results, 
£ t 2 0Ct is negative and finite and remains so in the limit 
in which one or more dimensions of the rectangle vanish 
and up to (d — 1) pairs of hyper-planes coincide. Note 
that £^ cct (£k 0) = ^£ r 2 oct . ■ When any dimension of 
the rectangle becomes large £j? cct vanishes. As mentioned 
previously this analyticity in the shape and dimensions of 
the objects is expected in the world-line description and 
is one of the more interesting characteristics of the irre- 
ducible iV-body Casimir energies defined by eq. © and 
eq. (JSJ) - One might ask whether eq. (|19l) can be given a 
meaning that does not involve divergent irreducible con- 
tributions involving less than 2 d -bodies. In fact £ r 2 0Ct de- 
termines higher order derivatives of the vacuum energy 
in much the same manner as the original Casimir energy 
determines the force between two plates: 



dp, dp 



■ dtjErert,. — Of. dp. 



(20) 



since all irreducible contributions with less than 2 d bodies 
do not depend on all dimensions of the rectangle. For 
d = 2 pairs of parallel plates eq. (|20|) implies that the 
irreducible four-body contribution to the vacuum energy 
fully describes certain stability derivatives of the vacuum 
energy. Other irreducible finite iV-body Casimir energies 
also suffice to describe, generally more complicated, higher 
order derivatives of the vacuum energy. Finite irreducible 
iV-body Casimir energies do not depend sensitively on the 
(quantum) description of the intersecting objects at high 
scattering energy and thus are reliably estimated by low- 
energy effective models that describe the interaction with 
the objects in terms of potentials or boundary conditions. 

An estimate of the magnitude of irreducible three-body 
Casimir forces is provided by the Casimir- Polder force on a 
polarizable atom due to a bilayer. Using the results of [5T], 
fig. [5] compares the irreducible electromagnetic three-body 
Casimir force on an atom with the sum of irreducible two- 
body Casimir forces for a Si/Si02 bilayer. It is well known 
that Casimir forces are not additive and the three-body 
correction in this example is not negligible. It diminishes 
the overall attractive force by almost 15% for atoms that 
are about 10 layer-thicknesses from the bilayer. 

Irreducible Casimir Energies in the Multiple 
Scattering Expansion. — To address more compli- 
cated geometries in the electromagnetic case, a represen- 
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Fig. 2: Ratio of the irreducible three-body Casimir force to the 
sum of irreducible two-body Casimir forces on an atom near a 
Si/Si02 bilayer. The distance of the atom from the surface is 
measured in units of the thickness, b, of the SiC>2 layer. 

tation of the irreducible Casimir energies in terms of one- 
body T-matrices is required [T2]. For three bodies it may 
be obtained using the generating functional approach of 
[2"9"] . The irreducible three-body Casimir energy 

£(3) 

ex- 
pressed in terms of the free-, 1-, 2- and 3-body Green's 
functions in the notation of 113) is, 



?(3) 



^Tr(ln G 123 - In G 12 - In G 23 - In G 13 
+ In Gi + In G 2 + In G 3 - In G$) (21) 
gTrOnGiG^Gaa - InGiG^Ga - In GxG^G z ) , 



where G a = G%G a is the Green's function for the domain 
2) Q . The trace is over space and time, with r here denoting 
the temporal extent. Using G" 1 = G^ 1 + GJ 1 — 1 and 



^123 = G i 1 + G 23 ~ 1 with Gi = l-Ti, the irreducible 
three-body Casimir energy of eq. (|21j) in terms of one-body 
scattering matrices Tj finally is, 

£(3) = ^TrHU-fi^-Gaa^+lnM-lnlXia] 
= ^Tr In [l - A 12 f x (f^fa - G 2 f 3 X 23 f 2 

— G 3 T 2 A 32 T 3 ^ X13 . (22) 

Here Tj = T; L G% = (1 — G;), with G% the Green's function 
for the domain ®q with no objects inserted. The operators 
Xij satisfy the integral equation, 



Xij(1 — fifj) = 1 



(23) 



The expression in eq. (|22j) differs from that given in [12] 
only in that all two-body interactions have been sub- 
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tracted. As our previous considerations show, the irre- 
ducible three-body Casimir energy given ineq. (j2"2")l is con- 
tinuous in the position of the three objects and remains 
finite when two overlap pairwise even as the corresponding 
(irreducible) two-body contribution to the vacuum energy 
diverges. Every term in eq. (|22[) requires scattering off all 
three objects and is individually finite. We explicitly cal- 
culated [28 the three-body correction of eq. (1221) to the 
Casimir energy of three semi-transparent parallel plates 
[28] . The result is symmetric under the exchange symme- 
try and remains finite when any two of the three plates 
coincide. Writing irreducible iV-body Casimir energies in 
terms of scattering matrices unambiguously defines them 
for any local field theory and in particular for the electro- 
magnetic case. 

* * * 
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